Violation of the Minimum H-H Separation "Rule" for Metal Hydrides 
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Using gradient-corrected, all-electron, full-potential, density-functional calculations, including 
structural relaxations, it is found that the metal hydrides .RTInHi.333 (R — La, Ce, Pr, or Nd; 
T = Ni, Pd, or Pt) possess unusually short H— H separations. The most extreme value (1.454 A) 
ever obtained for metal hydrides occurs for LaPtInHi.333. This finding violates the empirical rule for 
metal hydrides, which states that the minimum H— H separation is 2 A. Electronic structure, charge 
density, charge transfer, and electron localization function analyses on .RTInHi.333 show dominant 
metallic bonding with a non- negligible ionic component between T and H, the H— H interaction 
being weakly metallic. The paired, localized, and bosonic nature of the electron distribution at the 
H site are polarized towards La and In which reduces the repulsive interaction between negatively 
charged H atoms. This could explain the unusually short H— H separation in these materials. Also, 
R—R interactions contribute to shielding of the repulsive interactions between the H atoms. 

PACS numbers: 71., 81.05.Je, 71.15.Nc, 71.20.-b 



The most attractive aspect of metal hydrides from 
a technological point of view is their potential use as 
energy-storing materials. Hydrogen as energy carrier and 
other visions of the "hydrogen society" are especially at- 
tractive from an environmental point of view, but since 
hydrogen is a low-density gas at STP, the storage of the 
large quantities required for most applications is a chal- 
lenge. High-pressure-compressed-gas storage is energy 
intensive if high volume efficiency is desired, liquid or 
solid hydrogen storage even more so, and all involve cer- 
tain hazards. Storage of hydrogen in the form of solid 
metal hydrides, from which it can readily be recovered 
by heating, is safe and volume efficient. 

The amount of hydrogen per volume unit in metal hy- 
drides is very high; in some cases higher than in liquid 
or solid hydrogen, e.g., VH2 stores more than twice the 
amount of hydrogen than solid H2 at 4.2 K. It is unfor- 
tunate, however, that most metal hydrides are heavy in 
relation to the amount of hydrogen they contain. FcTiH2 
and LaNi5H7, e.g., only contain 1.9 and 1.6wt.% hydro- 
gen, respectively. Therefore, efforts in hydride research 
have over the past 25—30 years been concentrated on de- 
signing new, or modifying known, intermetallic hydrides 
to increase the storing capacity and simultaneously ad- 
justing their properties to make them capable of deliver- 
ing hydrogen at useful pressures (>0.1 MPa) and accept- 
able temperatures (< 425 K) These aspects are par- 
ticularly important for most mobile applications where 
hydrogen would be used directly in combustion engines 
or indirectly via fuel cells. It has proven difficult to ex- 
ceed 2 wt.% of stored hydrogen and it remains a challenge 
to increase this figure if metal hydrides are to become a 
viable source for the transportation sector. 

The search for efficient hydrogen-storage metal hy- 
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drides g] has to some extent been hampered by the 
mental barriers which empirical rules have put on the 
thinking. For example, the interstitial hole size where 
hydrogen is expected to occupy should be > 0.40 A. Swi- 
tendick || observed from a compilation of experimental 
structure data that the minimum H— H separation in or- 
dered metal hydrides is >2A ("the 2- A rule"). This em- 
pirical pattern is later Q supported by band-structure 
calculations which ascribe the effect to repulsive interac- 
tion generated by the partially charged hydrogen atoms. 
A practical consequence of this repulsive H— H interac- 
tion in metal hydrides is that it puts a limit to the amount 
of hydrogen which can be accommodated within a given 
structural framework. So, if H— H separations less than 
2 A would be possible this could open for new efforts to 
identify potential intcrmctallics for higher hydrogen stor- 
ing capacity. However, there are indeed metal hydrides 
which do violate "the 2-A rule" and we have here identi- 
fied the origin for such behavior. 

i?NiIn (R = La, Ce, Pr, and Nd) crystallizes in the 
ZrNiAl-type structure (space group PQ2m) and can for- 
mally be considered as a layered arrangement with a re- 
peated stacking of two different planar nets of composi- 
tion i?3Ni2 and Niln3 along [001] of the hexagonal unit 
cell. If the hydrides of these materials obey the hole-size 
demand and "the 2-A rule" one would expect H to occupy 
the interstitial 2d site within .R3M2 trigonal bipyramid. 
However, proton magnetic resonance (PMR) studies sug- 
gest |5|, |(| that H occupies both the Ah and 6i sites or 
either of them with H— H distances in the range 1.5—1.8 
A. Recent powder X-ray and neutron diffraction stud- 
ies [0 on i?NiInDi.333_a; (ideally i?3Ni 3 In 3 D4) show that 
deuterium occupies the Ah site located on three-fold axis 
of i? 3 Ni tetrahedra that share a common face to form 
trigonal bipyramid (Fig.|l|). This configuration gives rise 
to extraordinary short H— H separations of around 1.6 
A 0. As the diffraction techniques generally determine 
the average structure, neglect of partial H-site occupan- 
cies and local lattice distortions may lead one to conclude 
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FIG. 1: The i?NiInHi. 3 33-type crystal structure. The i? 3 Ni 2 
bipyramid is emphasized with thicker lines. 



with shorter H— H separations than actually present in 
the real structure [||. Hence, it is of interest to perform 
structural optimization theoretically. 

The present full-potential LMTO j| calculations are 
all electron, and no shape approximation to the charge 
density or potential has been used. The basis functions, 
charge density and potential were expanded in spheri- 
cal harmonic series inside the muffin-tin spheres and in 
Fourier series in the interstitial regions. The ratio of 
interstitial to unit cell volume is around 0.42. The calcu- 
lations are based on the generalized-gradient-corrected- 
dcnsity- functional theory as proposed by Perdew et 
al. flCfl Spin-orbit terms are included directly in the 
Hamiltonian matrix elements for the part inside the 
muffin-tin spheres. The basis set contained semicore hp 
and valence 5d, 6s, 6p, and 4/ states for La (for Ce the 
4/ electrons are treated alternatively as valence and lo- 
calized core electrons, whereas the Nd-4/ and Pr-4/ elec- 
trons are treated as localized electrons using open core 
approximation), 3s, 2>p, 4s, Ap, and 3d for Ni, 4s, Ap, 5s, 
5p, and Ad for Pd, 5s, 5p, 6s, 6p, and 5d for Pt, 5s, 5p, 
and 5d for In and Is, 2p, and 3d states for H. All orbitals 
were contained in the same energy panel. To ensure well- 
converged wave functions a so-called multi basis was in- 
cluded, implying the use of different Hankel or Neuman 
functions each attaching to its own radial functions. This 
allows reliable description of the higher-lying unoccupied 
states and low-lying semicore states. The self consistency 
was obtained with 105 k points in the irreducible part of 
the Brillouin zone. Test calculations were made for 210 
k points to check for convergence, but the optimized c/a 
ratio for 105 and 210 k points for LaNiInHi.333 is essen- 
tially the same. To gauge the bond strength and nature 
of bonding we have used crystal orbital Hamiltonian pop- 
ulation (COHP) and electron localization function (ELF) 
analyzes, as is implemented in TBLMTO-47 p~T[ | . 

All calculations relate to ideal and fully saturated hy- 
drides with composition i?3T3in3H4 (flTInHi.333; R = 
La, Ce, Pr, or Nd, T = Ni, Pd, or Pt). For (trivalent) R 
the 4/ electrons were treated as core electrons (except for 



La; 4/). As Ce-4/ electrons are known to take different 
valence states in intermetallic compounds, different pos- 
sibilities Q for the valence states of Ce were considered 
in both the hydrides and intermetallic phases during the 
structural optimization. For this optimization all atom 
positions were relaxed by force minimization and equi- 
librium c/a and volume were obtained by total energy 
minimization. Optimized structural parameters for the 
Ce compounds are in good agreement with experimental 
values (only) when Ce atoms are assumed to be in the 
trivalent state. The calculated equilibrium lattice pa- 
rameters and the changes between the intermetallic and 
corresponding hydride phases are given for selected com- 
pounds along with experimental parameters in Table Q. 

In general, the calculated lattice parameters are in 
good agreement with the experimental values, and the 
small differences found may partly be attributed to hy- 
drogen non-stoichiometry (around 10%) in the experi- 
mental studies. The hydrogen-induced lattice expansion 
is strongly anisotropic (Table ||); a huge expansion along 
[001] (Ac/c = 14—20%) and a smaller intralayer con- 
traction (— Aa/a = 0—5.8%). The calculated cohesive 
energy and heat of formation for the hydrides are larger 
than for the corresponding intermetallic phases indicat- 
ing that it might be possible to synthesize all these hy- 
drides. The electronic structure studies show that all con- 
sidered phases are in the metallic state consistent with 
experimental findings. The calculated R— H, T— H, and 
H— H distances are given in Table H along with experi- 
mentally available values. An interesting observation is 
that all i?TInHi.333 materials have unusually short H— H 
distances. Two explanations have been proposed. Pair- 
ing of the hydrogen atoms (either by molecular H2-like 
bonding or by bonding mediated by the intermediate T 
atom) has been advanced to explain the anomalous PMR 



spectrum of CeNilnHi.o 
focuses on the significant 



13|. The second explanation 
y shorter La— La distance in 
LaNiInHi.333 than in closely related phases Q, whereby 
the La atoms (generally R) may act as a shielding that 
compensates the repulsive H— H interaction. 

In order to evaluate these possibilities we have calcu- 
lated the total energy for (hypothetical) LaPtlnHi 333 as 
a function of H— H separation according to three different 
scenarios: (1) Keeping La, Pt, and In fixed in their equi- 
librium positions, (2) moving La 0.08 A out of the equi- 
librium position toward H, and (3) moving La 0.08 A out 
of the equilibrium position away from H. The obtained 
results are illustrated in Fig.|[ When La, Pt, and In are 
in their optimized equilibrium position, the equilibrium 
H— H separation is 1.454 A. This scenario corresponds to 
a lower total energy than the two alternatives. For sce- 
nario 3 we obtain a shorter H— H separation (1.438 A) 
than for the ground state configuration, and for scenario 
2 a correspondingly larger separation (1.462 A). 

As the total energy curves increase steadily on reduc- 
tion of the H— H separations, the possibility of stabiliza- 
tion of hydrogen in the form of molecular H2-like units 
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TABLE I: Calculated lattice parameters (a and c in A) and c/a for LaTInHi.333 and relative variation in unit cell dimensions 
(in %) consequent on hydrogenation from RTIn to -RTInHi.333. 



Compound 


a 


c 


c/a 


Aa/a 


Ac/ 


c 


AV/V 


The. 


Exp. 


The. 


Exp. 


The. 


Exp. 


The. 


Exp. 


The. 


Exp. 


The. 


Exp. 


LaNilnH 1.333 


7.2603 


7.3810 


4.5522 


4.6489 


0.6270 


0.6399 


-3.969 


-2.76 


14.02 




14.8 


5.14 


8.54 


LaPdlnH 1.333 


7.3501 




4.8112 




0.6546 




-5.42 




16.64 






4.33 




LaPtInHi.333 


7.7274 




4.6903 




0.6070 




-0.04 




13.98 






14.00 




CeNiInHi.333 


7.4536 


7.2921 


4.4871 


4.6238 


0.6020 


0.6341 


-1.68 


-3.21 


12.72 




16.3 


8.97 


8.98 


PrNiInHi.333 


7.3783 


7.260 


4.4726 


4.560 


0.6062 


0.6281 


-2.85 


-3.73 


13.93 




15.4 


7.52 


7.01 


NdNiInHi.333 


7.2408 


7.2255 


4.5560 


4.5752 


0.6292 


0.6332 


-3.72 


-3.92 


16.75 




16.5 


7.60 


7.53 



TABLE II: Calculated interatomic distances (in A) for 
firinHi.333. 





R- 


H 


T 


-H 


H 


-H 


Compound 


The. 


Exp. 


The. 


Exp. 


The. 


Exp. 


LaNilnH 1.333 


2.379 


2.406 


1.489 


1.506 


1.573 


1.635 


LaPdInHi.333 


2.373 




1.644 




1.523 




LaPtInHi.333 


2.475 




1.618 




1.454 




CeNiInHi.333 


2.427 


2.371 


1.457 


1.508 


1.572 


1.606 


PrNilnH 1.333 


2.387 




1.492 




1.487 




NdNiInHi.333 


2.350 


2.350 


1.493 


1.506 


1.492 


1.562 




1.4 1.5 
H-H distance ( A) 

FIG. 2: Total energy versus H— H distance in LaPtInHi.333. 
(a) All atoms except H are fixed at their equilibrium posi- 
tions, (b) La atoms are moved 0.08 A out of their equilibrium 
position toward H. (c) La atoms are moved 0.08 A out of their 
equilibrium position away from H. 



seems completely ruled out. The total energy increases 
drastically also for increased H— H separation beyond 
the equilibrium value. This is due to increasing repul- 
sive T— H interaction and decreasing attractive H— H in- 
teraction. The considerable changes in the equilibrium 
H—H distance on R displacement indicate that R in the 
R3T2 trigonal bipyramidal configuration (Fig.Q) acts as a 
shielding that to some extent compensates repulsive H—H 
interactions. When the R— H separation is reduced the 



H atoms are allowed to approach each other more closely. 

Owing to charge transfer from metal to hydrogen the 
repulsive H—H interaction in metal hydrides are generally 
larger than that within the H2 molecule and this may be 
the physical basis for "the 2-A rule" . Another example of 
violation of "the 2-A rule" is found for Th 2 AlH 4 |j, but 
here the minimum H—H separation (1.945 A) within the 
TI15 trigonal bipyramids is much larger than the bond 
distance in the H2 molecule and the H—H separations in 
i?rinHi.33 3 . In order to quantify the bonding interaction 
between the constituents in the i?TInH! 333 series the in- 
tegrated crystal orbital Hamilton population (ICOHP) 
were calculated. For example, the ICOHP values up to 
E F for LaNiInHi. 333 are -3.44, -0.14, -0.72, -0.85, 
-0.86, -1.21, and -0.61 eV for Ni(2c)-H, H-H, La-H, 
Ni(2c)-In, La-Ni(lb), Ni(lb)-In, and Ni(2c)-La, re- 
spectively. This indicates that the strongest bonds are 
those between Ni(2c) and H. Another important obser- 
vation is that the bonding interaction between the hydro- 
gens is small, which further confirms that the short H—H 
separation in these materials arc not rooted in hydrogen 
pairing or formation of H 2 -like molecular units. 

In order to substantiate this observation further we 
have calculated the valence-charge-density distribution 
in (100) of LaNiInHi.333 (Fig. [§>,). From this figure, it is 
apparent that Ni and H form an N1H2 molecule-like struc- 
tural subunit. Moreover, Fig.||a demonstrates that there 
is no substantial charge density distributed between the 
H atoms. In order to depict the role of charge transfer we 
have displayed the charge transfer (the difference in the 
electron density of the compound and that of constituent 
atoms superimposed on the lattice grid) for LaNiInHi.333 
within (100) in Fig. |b. From Fig. Elb it is clear that elec- 
trons are transfered from La, In, and Ni to the H site. So, 
there is considerable ionic bonding component between 
H and the metallic host lattice. The transferred electrons 
from the metallic host lattice to the H2-like subunit of the 
structure enter the antibonding a* levels and gives rise to 
repulsive interaction. This repulsive interaction between 
the negatively charged H atoms could explain why the 
H-H separation in these materials are larger than that 
in the H 2 molecule. If there was strong covalent bond- 
ing between Ni and H one should expect a significant 
(positive) value of charge transfer distribution between 
these atoms (contributed by both atoms). The absence 
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FIG. 3: (a) Total charge density, (b) Charge transfer and (c) Electron localization function plot for LaNiInHi.333 in the (100) 
plane. The origin is shifted to (1/3, 0, 0) and the charge densities are in e/a.u. 3 . 



of such a feature rules out this possibility. The ELF is 
an informative tool to distinguish different bonding in- 
teractions in solids pl| and ELF for LaNiInHi.333 in 
(100) is given in Fig.gc. The large value of ELF at 
the H site indicates strongly paired electrons with lo- 
cal bosonic character. Another manifestation of covalent 
bonding between Ni and H should have been paired elec- 
tron distribution between these atoms. The negligibly 
small ELF between Ni and H indicates that the proba- 
bility of finding parallel spin electrons close together is 
rather high (correspondingly small for antiparallel spin 
pairs) in this region confirm metallic bonding consistent 
with the result obtained from charge transfer analysis 
and the detailed analysis show that delocalized metallic 
Ni(2c)-d electrons are distributed in this region. Even 
though the charge distribution between Ni and H looks 
like a typical covalent bonding the charge transfer and 
ELF analyzes clearly show that the electron distribution 
between Ni and H are having parallel spin alignment and 
purely from Ni site. Hence, chemical bonding between Ni 
and H is dominated by metallic components with consid- 
erable ionic weft. The partial density of state analysis 
also show that the H-s states are well separated from 
the Ni-d states in the whole valence band indicates the 



presence of ionic bonding between Ni and H. Due to the 
repulsive interaction between the negatively charged H 
electrons, the ELF contours are not spherically shaped 
but polarized towards La and In. The localized nature 
of the electrons at the H site and their polarization to- 
wards La and In reduce significantly the H—H repulsive 
interaction and this can explain the unusually short H—H 
separation in this compound. The ELF between the H 
atom takes a significant value of 0.35. Considering the 
small charge density, this indicates a weak metallic type 
of interaction between the hydrogen atoms. 

i?TInHi.333 constitutes a series with much shorter 
H—H separations than other known metal hydrides. We 
have shown that the short distances between the H atoms 
in such metal hydrides are governed primarily by the 
polarization of negative charges on H towards the elec- 
tropositive La and In. We believe that this conclusion 
is of more general validity, and may be utilized to search 
for other metal hydrides of potential interest as hydrogen 
storage materials. 

The authors are grateful to the Research Council of 
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the Norwegian supercomputer facilities. 



[1] A.J. Maeland in Recent Advances in Hydride Chemistry, 
edited by R. Poli (North-Holland, Amsterdam, in press). 

[2] L. Schlapbach, F. Meli and A. Ziittel in Intermetalhc 
Compounds: Practice, vol. 2, edited by J.H. Westbrook 



and R.L. Fleischer (Wiley, New York, 1994) pp. 475. 
[3] A.C. Switendick, Z. Phys. Chem. 117, 89 (1979). 
[4] B.K. Rao and P. Jena, Phys. Rev. B 31, 6726 (1985). 
[5] K. Ghoshray, B. Bandyopadhyay, M. Sen, A. Ghoshray, 



■5 



and N. Chatterjee, Phys. Rev. B 47, 8277 (1993). 

[6] M. Sen, A. Ghoshray, K. Ghoshray, S. Sil, and N. Chat- 
terjee, Phys. Rev. B 53, 14345 (1996). 

[7] V.A. Yartys, R.V. Denys, B.C. Hauback, H. Fjellvag, I.I. 
Bulyk, A.B. Riabov, and Ya. M. Kalychak, J. Alloys and 
Comp. 330-332, 132 (2002). 

[8] K. Yvon and P. Fischer in Hydrogen in Intermetallics, 
edited by L. Schlapbach (Springer, Berlin, 1988). 

[9] J. M. Wills, O. Eriksson, M. Alouani, and D. L. Price 
in Electronic structure and Physical Properties of Solids, 
edited by H. Dreysse (Springer, Berlin, 2000) p. 148. 
[10] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 



[11] G. Krier, O. Jepsen, A. Burkhardt, and O. K. Ander- 
sen, Tight Binding LMTO-ASA Program Version 4.7 
(Stuttgart, Germany, 2000). 

[12] A. Delin and B. Johansson, J. Magn. Magn. Mater. 177- 
181, 373 (1998). 

[13] K. Ghoshray, B. Bandyopadhyay, M. Sen, A. Ghoshray, 
and N. Chatterjee, Phys. Rev. B 47, 8277 (1993). 

[14] P. Vajeeston, R. Vidya, P. Ravindran, H. Fjellvag, A. 
Kjekshus and A. Skjeltorp, Phys. Rev. B 65, 075101 
(2002). 

[15] A. Savin, R. Nesper, S. Wengert, T. Fassler, Angew. 
Chem. Int. Ed. Engl, 36, 1809 (1997). 



